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Abstract: We define and study a set of operators that compute statistical properties 
of the Brownian Loop Soup, a conformally invariant gas of random Brownian loops 
(Brownian paths constrained to begin and end at the same point) in two dimensions. 
We prove that the correlation functions of these operators have many of the properties 
of conformal primaries in a conformal held theory, and compute their conformal dimen¬ 
sion. The dimensions are real and positive, but have the novel feature that they vary 
continuously as a periodic function of a real parameter. We comment on the relation 
of the Brownian Loop Soup to the free held, and use this relation to establish that the 
central charge of the Loop Soup is twice its intensity. 
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1 Introduction 

1.1 The Brownian Loop Soup 

Take a handful of loops of various sizes and sprinkle them onto a flat surface. The 
position where each loop lands is uniformly random, independent of any loops already 
in place. Each loop is Brownian - a Brownian motion constrained to begin and end at 
the same “root” point, but otherwise with no restriction - and characterized by a “time” 
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Figure 1: Two Brownian loops, each of time length t = 1. 

length t that is linearly related to its average area (c/. Fig. 1). The distribntion in t is 
~ dt/t^, so that there are many more small loops than large, and is chosen to ensure 
invariance under scale transformations. The overall density of loops is characterized by 
a single parameter: the “intensity” A > 0. This random ensemble of loops is called the 
Brownian Loop Soup (BLS) and was introduced in [1]. 

More precisely, the BLS is a Poissonian random collection of loops in a planar 
domain D with intensity measure where A > 0 is a constant and is the 

restriction to D of the Brownian loop measure 



( 1 . 1 ) 


where A denotes area and is the complex Brownian bridge measure with starting 
point z and duration t. We note that the Brownian loop measure should be interpreted 
as a measure on “unrooted” loops, that is, loops without a specihed “root” point. 


(Formally, unrooted loops are equivalence classes of rooted loops—the interested reader 
is referred to [1] for the details.) For ease of notation, the /i^°°^-measure of a set {...} 
will be denoted ..) = ..}). 


The BLS turns out to be not just scale invariant, but fully conformally invari¬ 
ant. For sufficiently low intensities A, the intersecting loops form clusters whose outer 
boundaries are distributed like Conformal Loop Ensembles (CLEs)[2]. CLEs are the 
unique ensembles of planar, non-crossing and non-self-crossing loops satisfying a natural 
conformal restriction property that is conjecturally satished by the continuum scaling 
limits of interfaces in two-dimensional models from statistical physics. The loops of a 
CLEk are forms of SLE^ (the Schramm-Loewner Evolution with parameter k [3]). The 
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Figure 2: The Brownian loop soup for intensities A = 2.5 (left panel) and A = .5 (right 
panel). The plots were made with both a “UV” cutoff on short timelengths and an “IR” 
cutoff on long timelengths. The true BLS is a scale-invariant fractal that covers every point 
with probability 1. 

CLEs generated by the BLS correspond to values of k between 8/3 and 4. For example, 
the collection of outermost interfaces in a planar critical Ising model in a hnite domain 
with plus boundary condition is conjectured to converge to CLE 3 in the scaling limit. 

In this paper we will dehne and compute certain statistical correlation functions 
that characterize aspects of the BLS distribution. We will focus on two types of in¬ 
formation: the number of distinct loops that surround (or “cover” if one thinks of the 
loops as being hlled in) a given point or set of points, and the net number of windings 
of all the loops around a given point or points (see Fig. 3). In both cases, we hnd 
results consistent with the correlation functions of primary operators in a conformal 
held theory. 

1.2 Motivation 

We have several motivations for this work. In [4] one of us considered a similar model, 
where instead of Brownian loops one sprinkles disks. The 2- and 3-point correlation 
functions of certain operators in that model behave like those of a conformal held theory 
(GET), and with a novel set of conformal dimensions. This is of interest because the 
model of [4] was derived in [5] as an approximation to the asymptotic distribution of 
bubble nucleations in theories of eternal inhation, a theory for which there is some 
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Figure 3: A stylized Brownian loop. The numbers indicate the winding numbers of the 
loop that contribute additively to N^,, while the green shaded region is the interior of the 
loop (the set of points disconnected from infinity by the loop) that contributes ±1 (where the 
sign is a Boolean variable assigned randomly to each loop) to the layering number Ni. 

reason to believe a CFT dual may exist [6-8]. However, [4] computed the 4-point 
function exactly, and it suffered from a deficiency: it was not smooth as the position 
of the fourth point crossed the circle connecting the other three. 

The origin of the non-analyticity in [4] is probably the fact that although the disk 
distribution is invariant under global conformal transformations, it is not locally con¬ 
formally invariant (since disks do not map to disks). By contrast the BLS distribution 
is fully conformally invariant. Therefore we expect the analogous correlation functions 
to be better behaved, perhaps dehning a healthy CFT that could be related to the 
physics of de Sitter spacetime and eternal inflation. 

Another motivation for considering CFTs related to the BLS is the relation between 
the BLS and SLE [3], which in turn is related to a large class of conformally invariant 
models ranging from percolation [9-12] to the Ising model [13]. 

If in fact the BLS correlation functions we study arise from or dehne a CFT, it 
does not seem to be one that is currently known, and it has several interesting and 
novel features. As in [4], the conformal dimensions of the primary operators are real 
and positive, but vary continuously and as a periodic function of a real parameter (3. 
As we will see this periodicity arises because the operators are of the form where 
N is integer valued. 

1.3 Correlation functions in the Brownian Loop Soup 

As mentioned above, we will study the correlation functions {i.e. expectation values) 
of two distinct types of operators in the BLS. 

Layering number: The hrst type is closely related to the operators considered in [4] 
(which we refer to as the “disk model”). For each loop, dehne the interior as the set 


4 




of points inside the outermost edge of the loop (including any isolated “islands” that 
might appear inside due to self-intersections; cf. Fig. 3). Points in this set are “covered” 
by the loop. For a point z, consider an operator N{z) that counts the number of distinct 
loops that cover the point z, so that N{z) is the number of “layers” at 2 ;. 

One difficulty arises immediately. Because the BLS is conformally (and therefore 
scale-) invariant, any given point of the plane is covered by inhnitely many loops with 
probability one. Since N{-) > 0, this means that, for any fixed z, {N{z)) diverges with 
probability one. Precisely the same difficulty arises in the disk model, where it was 
dealt with by adding another, identical and independent copy of the distribution, and 
then counting the difference in the number of disks of each type that cover z.^ We will 
follow a very similar procedure here, assigning to each loop a random Boolean value 
and then defining the layering operator Ni{z) = N+{z) — N_{z). 

Winding number: The other operator we will discuss is N^{z), which counts the 
total number of windings of all loops around a point z (Fig. 3). This makes use of the 
fact that Brownian loops have an orientation (they grow in a particular direction as 
the time t increases). Since the winding number can be positive or negative, it is not 
necessary to include another copy of the distribution or compute a difference of two 
values. 

A natural physical interpretation of winding in the BLS is as follows. If each loop 
represents the conhguration of a string at some instant of time in a 2-1-1 dimensional 
spacetime, winding number counts the number of units of the flux the string is charged 
under (see e.g. [15]). 

In both cases we will focus on exponentials of these operators times imaginary 
coefficients (correlators of the number operators themselves are afflicted by logarithmic 
divergences, like massless fields in two dimensions). Due to the similarity to free-held 
vertex operators we will denote these as 

Vg{z) = 

where N can be either a layering or winding number operator.^ 

Our paper is written for a mixed audience of mathematicians and physicists. In 
most of the paper we present rigorous proofs of our results. Sections 5 and 6, where 
we perform “physics-style” calculations, are the exceptions. The appendix. Section 8, 

^In this disk model this can be thought of as counting the number of bubble transitions that have 
affected the point z, in a model where the field has a discrete shift symmetry and there is a simple 
rule for bubble collisions (namely the one discussed in [14]). 

special case of the winding model on a lattice was considered by Le Jan (see Section 6 of [16]). 
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is dedicated to two important lemmas about the Brownian loop measure (1.1) which 
are used several times in the rest of the paper. A reader uninterested in our methods 
may simply read Sec. 2 for a summary of results, and Sec. 7 for our conclusions. 


2 Summary and results 

Our main results relate to correlation functions (he. expectation values of products) of 
exponentials of the winding and layering operators in the BLS. Specihcally, we establish 
the following: 


• For both versions in hnite domains D, correlators of n G N exponential operators 






) 5,D 


( 2 . 1 ) 


exist as long as a short-time cutoff 5 > 0 on the loops is imposed, and 


lim 

(5-s.O 



nj=i 


0£,(2;i, ... ...,/?,,) = 0d(2;/3) 


exists and is hnite. Moreover, if D' is another hnite domain and f : D ^ D' is & 
conformal map such that z'^ = f{zi),... ,z'^ = f{zn), then 


= Y[\f{zj)\ (pniz-, f3), 
i=i 


where the A(/3) is dehned below. This is the behavior expected for a conformal 
primary operator. 


• For both versions in inhnite volume, correlators of n exponential operators 

= ( 2 . 2 ) 

vanish. However, in the case of the layering model, one can remove the short- 
time cutoh and still obtain a nontrivial limit by imposing the following “charge 
conservation” condition, satished mod 27r, 

'^Pj=27rk, keZ. (2.3) 

j 


• The correlators (2.2) of layering operators in the plane are hnite and non-zero 
when (2.3) is satished, so long as the loop soup is cut oh at short times S (no 
long-time cutoh is necessary). 
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In the case of 2 points, assuming (2.3), the 5 —0 limit of the renormalized 
layering operator correlators in the plane ( 2 . 2 ) can be explicitly computed up to 
an overall multiplicative constant. The result is 


0c(2:i, 2^2;/3i,/32) — C 2 


Zl - Z2 


A(/3i)+A(/32) 


where C 2 is a constant (Theorem 3.1). 

The (5 —)■ 0 limit of the renormalized 3-point function for the layering model in 
the plane, assuming (2.3), is 


(i)c{zi, Z2, Z2,] I3i, (32, (3^) = 


c. 


- 2 : 2 ! 


A(/3i)+A(/32)-A(/33) 


\Zl - ^ 31 


A(/3i)+A(/33)-A(/32) 


1^2 - 2 : 3 ! 


A(/32)+A(/33)-A(/3i) 


where C 3 is a constant (Theorem 4.5). 

The conformal dimensions A(/3) differ for the two types of operators. For the 
layering number. 


For the winding number. 


M/3) = ^(l-cos/3). 


A^(/3) = A/3(27r - I3)/M, 


where this formula is valid for 0 < /9 < 2ti, and A^(/3) is periodic under (3 —)■ 
(3 + 271. 


Open questions: Our results leave a number of questions to be answered. 

• We have not determined whether the correlators of the winding operators converge 
in the inhnite plane. The missing piece is the winding number analog of the fact 
that the BLS is thin [17] (see Sec. 3.4). It is plausible that these correlators do 
not in fact converge without an explicit long-distance regulator even when charge 
conservation (2.3) is satisfied. 

• We have not established that the n-point correlators in either model are analytic 
functions of the Zi, although we expect this is the case. 

• We expect that these correlators define some kind of conformal field theory. What 
conformal field theory is it? 
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• The disk model of [4] could be thought of as the late-time distribution of bubbles 
produced by a hrst-order phase transition in de Sitter spacetime. Is there an 
analogous physical interpretation of the BLS? 

3 Correlators of the layering and winding operators 

3.1 Correlators of layering operator 

As discussed in Sec. 1.3, Model 1 is dehned by randomly assigning a Boolean variable 
to each loop in the Brownian loop soup. Alternatively, one can think of this as two 
independent Brownian loop soups, each with a Poisson distribution with 

intensity measure where we take A+ = A_ = A/ 2 . (This follows from the 

fact that the collection of all loops from a BLS of intensity A+ and an independent one 
of intensity A_ is distributed like a BLS with intensity A+ -|- A_.) 

Denote by the number of loops 7 in the hrst (respectively, second) class 

such that the point 2 ; G C is separated from inhnity by the image of 7 in C. If 7 is 
the “£lled-in” loop 7 , then this condition becomes z & 'j, or z is covered by 7 . We are 
interested in the layering held Ni, with Ni{z) = N^{z) — N_{z). This is purely formal 
as both A^+(_)(z) are inhnite with probability one for any 2 ;. They are inhnite for two 
reasons: both because there are inhnitely many large loops surrounding 2 ; (infrared, or 
IR, divergence), and because there are inhnitely many small loops around z (ultraviolet, 
or UV, divergence). 

We will consider correlators of the exponential operator Vg = and show that 

there are choices of that remove the IR divergence and a normalization which removes 
the UV divergence. Specihcally, we are interested in the correlators V) 3 (zi,... ,Zn) = 
YYj=i ^i(u ) ~ and their moments 

(Ug( 2 'i, . . . , ^n)) = IE^(Vg(Zi, . . . , Z^)') 

where G C, /3 = (/3i,..., (3n) G and the expected value is taken with respect 
to the distribution Px^^ioop®Pii 2 -, where Px^^ioop is the Poisson distribution with intensity 
measure and P 1/2 is the Bernoulli distribution with parameter 1/2 (remember 

that each loop belongs to one of two classes with equal probability) or, equivalently, with 
respect to two independent copies of the Brownian Loop Soup with equal intensities 
A/2. 

As the held W has both IR and UV divergences, we get a meaningful dehnition 
by introducing cutohs which restrict the loops to have diameter^ within some 5 and 


^The diameter of a Brownian loop is the largest distance between two points on its outer boundary. 



R G 6 < R: let fisiii') = fl {7 : 5 < diam( 7 ) < R}) and consider the 

correlators 

{Vf,{z^, . . . , Zn))s,n ^ , 

where the expectation is as the above with replaced by 

3.2 The 1-point function in the layering model 

In this section we explicitly compnte the 1-point fnnction in the presence of IR and 
UV cntoffs. Replacing the area of a hlled Brownian loop of time length 1 with that of 
a disk of radins 1, the resnlt reprodnces the 1-point fnnction in the disk model [4]. 

Lemma 3.1. For all z E C, we have that 

/ R\~^d-cos^) 

w 

Proof. In this and other similar resnlts that follow, we hrst compnte {Vi 3 {z))s,r = 
Ea,(S,_r terms of probabilities snch as = fJ^ 5 ,R{'y ■ z E (see (3.1)). 

Snch calcnlations can be performed by nsing two independent copies of the BLS, 
with random distribntions of loops indicated by Mi and M 2 respectively; then writing 
N(_{z) = f ly^z(dMi — dM 2 ), where is the indicator fnnction of A] next, observing 
that the properties of the Poisson distribntion imply that Ea^^r _ 

g-/(i-e hnally, compnting the integral in this last expression. We now 

give a detailed proof based on the representation with randomly colored loops. 

Recall that with IR and UV cntoffs in place, the held Ni{z) can be realized as 
follows. Let 7 be a realization of loops, and let be a collection of independent 

Bernonlli symmetric random variables taking valnes in { — 1,1}. The qnantity 

Nt(z) = W s 

7Gr/,2:G7,<5<diam(7)</? 

is hnite Pa,ai*°°p almost snrely, since : z E 6 < diam( 7 ) < R} = fis,R{'y ■ z E 

7 } < 00 (see [18]). Now, 

(U/3(;^))5,K = EAAi? 

00 

= J2^x, 5,R Px,,s,ui^k), 

k=0 

where Ck = {v ■ \{l ^ V ■ ^ ^ < diam( 7 ) < R}| = k}. If X denotes a (±l)-valned 

symmetric random variable, 

Eaah = {E = (cos/?)^ 
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Therefore, for : z G ■j), we have that 



(3.1) 


_ g-Aa2,i,H(l-cos^) 


(3.2) 


Moreover, by Lemma 8.1 in Sec. 8 , 


Oiz,5,R = y iog T 
5 0 


1 R 


which implies 



A 

5 


(l-cos(/3)) 


{V^{z))s,R 


as claimed. □ 


3.3 The winding operator 

To dehne the second model, let N^i^z) denote the total winding number about the 


point z of all loops in a Brownian Loop Soup; as for the layering operators, this is 
a formal dehnition as Nw{z) might have divergences. Consider again the correlators 
Vi3{zi, ...,Zn) = and their moments {Vii{zi, ..., Zn)) = E{Vii{zi, ..., Zn)) 

where Zj & C, f3 = (/3i,... ,/9„) G M”, and the expected value is taken with respect 
to the BLS distribution. Denoting by P^^^ioop the Poisson distribution with intensity 
measure and restricting the loops to have diameter between some 5 and R G 


with 6 < R, we let fi’'sR{-) = fl {7 : 5 < diam( 7 ) < R}) and consider the 


correlators 



We now explicitly compute the 1-point function in the presence of IR and UV 
cutoffs for the winding model. 

Lemma 3.2. For all z E C, we have that 



(3.3) 


where the formula is valid for (3 G [0, 27r), and for (3 ^ [0,27r), in the right hand side, (3 
should he replaced by [(3 mod 27r). 
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Proof. For a point z and a loop 7 , let 0^{z) indicate the winding number of 7 around 
2 ;. Moreover, for fc G (N U {0})^ let 

= {v ■ |{7 ^ h • ^ ^ 7)^ ^ diam( 7 ) < R, \9^{z) \ = m}\ = km for all m G N}. 


If a loop 7 has \9.y{z) \ = m, then the winding number 0.y{z) is ±m with equal probability 
under P^^^ioop. Finally, using Lemma 8.2 from Sec. 8 , for m > 1 we have 

a,,5,R,m = k'tni'y : ^ ^ 7 , \9^iz)\ = m) 

= /i ^"°^(7 ■ z ^ ^,6 < diam( 7 ) < P, \9^{z) \ = m) 

1 R 
log-. 


TT^m^ ° (5 


For all fc G (N U {0})^ we have 


00 / •« \ Ic 

D I r \ _ n .-xo. 

^>',t^S.R\^k) ~ J_J_ ' • ^ 


z,S,R,m 


m=l 


K. 


as for different m’s the sets of loops with those winding numbers are disjoint. Hence 
with IR and UV cutoffs in place, denoting by E\^s,r the expectation with respect to 
the Poisson distribution ^loop with intensity measure we have, for all z, 


= E,_s,r 

= Y1 


fce(Nu{ 0 })'’ 


/ \ 




fce(Nu{o})f^ m=i 

-Aaz,5,H.m(l-cos(m/3)) 

m=l 


, /3(27r-/3) 


where the /3 on the right hand side of the last equality is to be interpreted modulo 27r. □ 


This coincides with the result (5.3) computed using physics path integral methods. 

3.4 The 2-point function in the layering model 

We now analyze the 2-point function when the IR cutoff is removed by the charge 
conservation condition (2.3). 
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Theorem 3.1. If jSi + /32 = 2k7r with /c G Z, there is a positive constant C 2 < oo such 
that, for all zi^ Z 2 , 


lim {Vfi,{zi)Vii^{z 2 ))^^ = C 2 

H^oo ’ 


1^1 - 2:21 


- ^ (2—cos /3i —cos /32) 


As a consequence, 


lim lim 




— 6*21^1 — Z 2 


- 4 (2—cos /3i- cos IS 2 ) 


l^OR-^oo (2-cos/3i-cos/32) 

Proof. Letting d= \zi — Zi\, for given (3i and [I 2 , and d> S, we have that 


where N 12 is the number of loops that cover both zi and Z 2 , and Ni {N 2 ) is the number 
of loops that cover zi but not Z 2 {z 2 but not zi, resp.). The two-point function factorizes 
because the sets of loops contributing to Nu, Ni and N 2 are disjoint; the S is replaced 
by d in the hrst factor in the second line because a loop covering both zi and Z 2 must 
have diameter at least d. 

As in the 1-point function calculation, we can write 

00 

^e*(/3i+/32)7Vi2^^^ _ ^ + / 32 ))”P;,^^*-p(iVi 2 =n) = e-^“^.«("n"2)(l-cos(/3i+/32))^ 

n=0 


where ad,R{zi,Z2) = fVdJil ■ zi,Z2 e 7 ). Similarly, if as,R{zi,^Z2) = ■ zi G 

7 , 2:2 ^ 7 ) and ttyij(“ 12 : 1 , 2 : 2 ) is correspondingly dehned, then 


_ -\as^ii(zi,^Z2){l-cosl3i) 

/S,R ^ 


and 

/„i02N2\ _ -Ao5_ij(-.zi,^2)(l-cos/32) 

/5,R ^ 


Combining the three terms we obtain 


> i(PiN({zi)+l32Ni{z2))\ 

\ /5,R 

_ g-Aad,fl(^i,22)(l-cos(/3i+/32)) g-Aa5_fl(zi,-.Z2)(l-cos/3i) g-AQ^_j{(-.zi,^2)(l-cos/32) 


It is easy to see that limj:j_,.oo ctid,R{zi, Z 2 ) = 00 . (This follows from the scale invariance 
of by considering an increasing - in size - sequence of disjoint, concentric annuli 
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around Zi and Z 2 that are scaled versions of each other.) Hence, in order to remove the 
IR cutoff, we must impose ( 2 . 3 ) and set /3i + /32 = 2 k 7 i, so that 1 — cos(/3i + /? 2 ) = 0. 
Assuming that /3i + ^2 = 2k7T, we are left with 

_ ^-\as, r(zi,^Z 2}{l-cos^i)-\as,Ri^zi,Z2)il-cos 132) 

To remove the infrared cutoff, we use the fact that the loop is thin: If zi ^ Z2, : 

e 7 , 2:2 ^ 7 , diam( 7 ) > 5) < 00 for any 5 > 0 (see [17], Lemma 4). By the obvious 
monotonicity of as^R^Zi, - 12 : 2 ) in R, this implies that 

lim as,R{zi, -. 2 : 2 ) = G 7 : 2:1 e 7 , 2:3 ^ 7 , diam( 7 ) > 5} = 05 ( 2 : 1 , -.^ 2 )- 

R —^00 

By scale, rotation and translation invariance of the Brownian loop measure 
05 ( 2 : 1 ,- 12 : 2 ) can only depend on the ratio x = d/ 6 , so we can introduce the notation 
a{x) = 05 ( 2 : 1 , - 1 ^ 2 )- The function o has the following properties, which are also imme¬ 
diate consequences of the scale, rotation and translation invariance of the Brownian 
loop measure. 

• a{x) = 05(0, -'z) for any z such that \z\ = d. 

• For cr > 1, if 5 < d, letting as,R{z) = az,5,R = t^&,R{.l : z e^), 

a{ax) = 05(0, -i(Tz) 

= 0^5(0, -^az) + oy^5(0, -^az) 

= a{x) + O5,tj5(0) = a{x) + Oi,cr( 0 ). ( 3 . 5 ) 


Now let 

G{x) = ^ ^^iWiNi{zi)+(32Ni{z2))'^^ . 

using (3.4) and the dehnition of the function o, we can write 

G[x) = g“^"(^)(2-cos/3i-cos/32) 

Then, for a > 1, 

G{ax) = 

Using Lemma 8.1, we have that 


ai,f7(0) = - logo. 
5 
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It then follows that, for cr > 1, 


G{ax) = 

For 0 < cr < 1, (3.5) implies 

a{ax) = a{x) — cri,i/cr(0). 


But since ^ 

tti,i/tT(0) = - - log cr, 

equation (3.6) is unchanged when 0 < cr < 1. 

The fact that (3.6) is valid for all cr > 0 immediately implies that 

G(a;) = 


(3.6) 


for some constant C 2 > 0. □ 

4 Conformal covariance of the n-point functions 

We now analyze the n-point functions for general n > 1 and their conformal invariance 
properties. In bounded domains D (Z C, we show, for both models, how to remove the 
UV cutoff (5 > 0 by dividing by , with the appropriate A/s. We also show 

that this procedure leads to conformally covariant functions of the domain D. The 
scaling with 6 originates from the fact that loops with diameter less than 6 can only 
wind around a single point in the limit 5 —)■ 0, and so for these small loops the n-point 
function reduces to the product of 1-point functions. 

In Section 4.3, we deal with the layering model in the full plane, C, and show that, 
together with the UV cutoff 5 > 0, we can also remove the IR cutoff R < 00 , provided we 
impose the condition ^ ^ttZ (c/. (2.3)). We refer to this condition as “charge 

conservation” because—apart from the periodicity—it is reminiscent of momentum or 
charge conservation for the vertex operators of the free boson. 

Just as for the two-point function, in the layering model the IR convergence (given 
“charge conservation”) is due to the hniteness of the total mass of the loops which 
cover some points but not others; this is basically the property that the soup of outer 
boundaries of a Brownian loop soup is thin in the language of Nacu and Werner [17]. 
We did not prove the analogous hniteness for the winding model, so in that case we 
are not able to prove that the IR cutoff can be removed. 
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4.1 The layering model in finite domains 


In the theorem below, we let (YYj=i^i3ji^j)'^ = ^x,s,d denote 

the expectation of the product 11^=1 with respect to a loop soup in D with 

intensity A > 0 containing only loops of diameter at least 5 > 0, that is, with respect to 
the distribution P^^ioop 0 P 1 / 2 , where P^^ioop is the Poisson distribution with intensity 

measure = /i^°^l{diam( 7 )> 5 } = h'''°^l{ 7 cD,diam( 7 )> 5 } and Pi /2 is the Bernoulli dis¬ 
tribution with parameter 1/2 (remember that each loop belongs to one of two classes 
with equal probability). 


Theorem 4.1. If n eN, D C C is bounded and (3 = {(3i,... ,(3n), then 


lllll -- 

(5-5>0 E/=i(1-cos/3P 


D 


(j)D{Zi, ...,Zn,l3) 


exists and is finite and real. Moreover, if D' is another bounded subset of C and 
f : D ^ D' is a eonformal map sueh that z'^ = f{zi ),..., 2 ;/ = f{zn), then 

n 

/3) = JJ \ f'{Zj)\-^^^-^°^^^^(j)D{zi, ...,Zn,f3). 

J=1 


The proof of the theorem will make use of the following lemma, where Bfiz) denotes 
the disc of radius 5 centered at 2 :, 7 denotes the complement of the unique unbounded 
component of C \ 7 , and where o(l) denotes a quantity “smaller than 0(1)”, i.e. that 
tends to zero as <5 —)■ 0 . 


Lemma 4.2. Let D,D' G C and let f : D ^ D' be a eonformal map. For n > 1, 
assume that zi,... ,Zn E D are distinet and that z'^ = f{zi ),..., 2 :/ = fizn), and let 
Sj = \ f'{zj) I for j = 1,..., n . Then we have that, for eaeh j = 1,... ,n, 

irfiTiG : e 7,7 f{Bs{zj)), 4 ^ 7 V/c ^ j) 

- ■■ z'^Efi,fi Bsjsiz'j), 4 ^ 7 V/c ^ j) = 0 ( 1 ) as 5 ^ 0. 

Proof. Let Bin{Zj) denote the largest (open) disc centered at Zj contained inside 
f{Bs{zj)) n Bs-siz'j), and Boufizj) denote the smallest disc centered at Zj containing 
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f{Bs{zj)) U A moment of thought reveals that, for 6 sufficiently small, 


e 7,7 ^ /(55 (zj)), 4 ^ 7 Vfc ^ j) 

- €■ Bs^siz'^), 4 ^ 7 Vfc ^ j)| 

= hD°^(7 : 4 e 7,7 f{Bs{zj)) n SjBs{z'j),^ C f{Bs{zj)) U SjBs{z'j), j) 

< : 4 e 7,7 5m(4),7 C BontiZj)) 


= clog 


diam(Eo„t(4)) 

diain{Bin ( 4 )) ’ 


(4.1) 


where c < cx) is a positive constant and the last equality follows from Proposition 3 of 
[18]. Note that, when D' = C, the quantities above involving are bounded because 
of the fact that the Brownian loop soup is thin [17]. 

Since / is analytic, for every w G dBs{zj), we have that 


\fM-Zj\ =Sj6 + 0{S^), 


which implies that 

diam(Eo„t(4)) _ 

/4o diam(Ei„(4.)) 

In view of (4.1), this concludes the proof of the lemma. □ 


Proof of Theorem 4.1. We hrst show that the limit is hnite. Using the notation of 
the previous section, we let rj denote a realization of loops and a collection of 

independent Bernoulli symmetric random variables taking values in { — 1,1}. Moreover, 
let [n] = { 1 ,..., u}, let /C denote the space of assignments of a nonnegative integer to 
each nonempty subset S of { 2 : 1 ,..., Zn}-, and for S C {zi, ..., Zn}, let Is C [u] be the 
set of indices such that /c G if and only if Zk G S. We have that 


I S,D 






keK. 


where Ck = {v ■ C {zi, ..., Zn}, S' 7 ^ 0, |{7 G 7 : G 7 Vj G Is, Zj ^7 Vj ^ 
Is}\ = k{S)}. With probability one with respect to have that for each 

j = l,...,n. 


Ne{z,) 


E E 





7:2jS7,diam(7)>(5 Sc{zi,...,Zn}:Zj£S 7:SC7,S'^C7'^,diam(7)>5 
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With the notation = J2r.sc^,s^cr4i^Mi)>s 
symmetric random variable, we have that 


, and letting X denote a (±l)-valne 


E 


\,S,R 



. 

Eaaj? E®(Eje/s|£^j 

Sc{^i,..., 2 „}, 5^0 

n (oos(Eft))‘‘^ 

SC{2l,...,2n},5/0 \ \jels / / 


Next, given S C {^i,...,^„} with 151 > 2, let a£)(5) = fioil : 5 C 7 , 5^^ C 7 '^) 
and ■ diam( 7 ) > S,Zj & Zk ^ ^ ^ j). Fnrthermore, let Is be 

the set of indices such that A; G Is if and only if Zk G S. Let m = minjj:j^j \zi — Zj\ f\ 
min* dist(zj, dD) and note that, when 6 < m, we can write 


k{S) 


=E n 


cos 


\j=l 


gD fceA; s'c{2i,.--,2„},|5|>i 


\k£ls 


{k{S))\ 


n(cos/?j 

i=i 


Mzj) XasM^j) 


n 

SC{21,...,2„},|S|>1 


(fc(^j))! 

-XaniS) I 1 - cos I /3fc 


\k&Is 


]^exp [-Xa5^D{zj){l - cos/dj)]. 
i=i 


(4.2) 


For every j = 1,..., n, using Lemma 8.1, we have that 


oc&,D{zj) = /i ^°^(7 : diam( 7 ) > 5, Zj e j,Zk ^ ^ j) 

= /i ^°^(7 : m > diam( 7 ) > 6, Zj G 7 ) 

+ Fd°^( 7 : diam( 7 ) > m, zj G 7 , 2 ;^ ^ 7 V/c 7 ^ j) 
1 7Tl 

= — + am,oi^j)- 

5 0 
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Therefore, we obtain 




5,D 


lim 

(5-s>0 Ej=i(l-cos(^j)) 


n 

SC{2l,...,2n},|S|>l 


-XaoiS) - cos 


1 E"=i(l-cos/3j)g- T ,]=1 Aam.oCzj)(1-cos/3j) 
= m~^ E"=i (1-cos/3j) 


exp 


-A 


S'C{21,...,2„},|S'| 

= (Pd{zi, ...,Zn;f3). 


-X'^am,D{Zj)il - cos(3j) 

J=1 

aD(^) |l-cos I j 

...,2„},|S'|>i V Vfce/s / . 


This concludes the hrst part of the proof. 

To prove the second part of the theorem, using (4.2), we write 


n 


n 



exp 


-A aB(^) 

5c{2i,...,2„},|5|>2 



]^exp [-Xas,D{zj){l - cos/^^)]. 
i=i 



(4.3) 


For each S C {zi,... ,Zn} with \S\ > 2 , an^S) is invariant under conformal transfor¬ 
mations, that is, if / : iA —?■ D' is a conformal map from D to another bounded domain 
D', and S' = {z(,...,4}, where z'^ = /(zi),...,^ = f{zn), then ao'iS') = aniS). 
Therefore, the hrst exponential term in (4.3) is also invariant under conformal trans¬ 
formations. This implies that, for <5 sufficiently small, 

n 

JJexp {-A [as^oizj) - as,D'{z'j)\ (1 - cos/3j)} . 
i=i 

Writing 



as,D{zj) = : diam(7) > 5 ,Zj G 7,7 C Bi,{zj)) 

+ hD °^(7 t Bs{zj),Zk ^ 7 Vfc ^ j) 

and noticing that /i ^°^(7 : diam(7) > 6,Zj G 7,7 C Bs{zj)) = /i^°^(7 : diam(7) > 
S, z'j G 7,7 C Bs{z'j)) (where we have assumed, without loss of generality, that 6 is so 
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small that Bs{zj) C D and Bs{Zj) C D'), we have that 

a5,D{zj) - as,D'iZj) = G 7,7 ^ Bs{zj), j) 

- : z' e 7,7 ^ Bsiz'j), j). 

To evaluate this difference, using conformal invariance, we write 

% e 7,7 ^ Bs{zj),Zk ^ 7 Vfc ^ j) 

= : 2 ;' e 7,7 f{Bs{zj)), 4 ^ 7 Vfc ^ j). 

Letting Sj = \f'{zj)\ and using (8.1) and Lemma 8.1 from the appendix, we can write 

O^S,D{Zj) — Oi5,D'{Zj) 

= : 4 ^ 7,7 fiBsizj)), 4 ^ 7 Vfc ^ j) - /iS°^(7 : 4 G 7,7 Bs{Zj), 4 ^ 7 V/c ^ j) 

= /^d°^(7 f{Bs{zj)), 4 ^ 7 Vfc ^ j) - /iB°^(7 : 4 G 7,7 B,^s{Zj), z'^^^^k^ j) 

- [/^d°^(7 : 4 G 7,7 Bs{z'j), 4 ^ 7 Vfc ^ j) - 4d°^( 7 : 4 G 7,7 ^s,5(4), 4 ^ 7 Vfc ^ j)] 

= /^S°^(7 ■ z'j^^n (t f{Bs{zj)), z'^i^Mk^ j) - /i^?^(7 : 4 ^ 7 ,7 44)’ 4 ^ 7 Vfc 4 j) 

- - logs,-. 

Using Lemma 4.2, we obtain 05 ,— C(s,D'{.Zj) = —| log \ f'{zj)\ + o(l) as h —)■ 0, 
which gives 


(14(^1, . . . , Zn))g jy 

( 14 (^ 1 , . . . , 4))5,D' 


= n If as h ^ 0. 

i=i 


Letting h —)■ 0 concludes the proof. □ 


4.2 The winding model in finite domains 


As above, let Nu,{z) denote the total number of windings of all loops of a given soup 
around z G C. We have the following theorem. 


Theorem 4.3. //n G N, H C C zs hounded and (3 = (/3i, ... ,/dn); then 

01 -in i^l) 


lim 

(5-s>0 


ft-) 


= 4d(U, ...,Zn,(3) 


exists and is finite and real. Moreover, if D' is another hounded subset of C and 
f : D ^ D' is a conformal map such that 4 = f{zi), ..., 4 = f{zn), then 

fiD'{z[, ..., 4; /3) = ]^ If (u)r^ ' ' 'ipnizi,..., z^, (3), 

j=i 

where, in the exponent, the fij’s are to he interpreted modulo 27i. 
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Proof. The proof is analogous to that of Theorem 4.1. Let a£){S; , ki^) := Hoin • 

6^{zi.) = ki. for each Zi- G S and C 7 '^)) and as,D{zj;k) := fioil ■ Zj G 'j,9^{zj) = 
k, Zk ^ ^ yk ^ j). With this notation we can write 


^ 0 nN'uj{Zn)^ 


5,D 


= exp 


-^E E E aoiS; fcii,..., fcij (1 - cos{ki^(3i^ + ... + 


1=2 SC{zi,...,Zn} ki. =-QO 

|S|=i 


n 

JJexp 

i=i 


-A as^D^Zj] k){l — cos{k/3j)) 


k=—oo 


For each S C {zi,...,Zn} with IS"! = / > 2, aoiS; ki^^,..., ki^) is invariant under 
conformal transformations; therefore, 


^gi/ 3 lWu,(^i) gi/JnW'iuizn)^ n C 00 

= E [o^s,D{^vk)-as,Diz^zk)]{l-cos{k(3,)) 

\ • • /5,D' j=l L k=-oo 

Proceeding as in the proof of Theorem 4.1, but using Lemma 8.2 instead of Lemma 8.1, 
gives 

a 5 ,D{zj; k) - as,D'{Zj-,k) = -c^log \f\zj)\ + o(l) as 5 -)■ 0 , 


where = 27 ^ for fc G Z \ {0} and Cq = 1/30. 

This, together with the observation, already used at the end of Section 3.3, that 
YlT=-oo ca:( 1 —cos(fc/?)) = (where, on the right hand side, /? should be interpreted 

modulo 27r), readily implies the statement of the theorem. □ 


4.3 The layering model in the plane 

Recall that (YYj=i^Pi{zj)\ denotes the expectation of the product ]/[? ^ 

\ / 6,R 

with respect to a loop soup in C with intensity A > 0 containing only loops 7 of 
diameter 0 < 5 < diam( 7 ) < R < 00 . 


Theorem 4.4. If n eN and /3 = (A, ■ ■ ■, /3n) with \(3\ = f^j ^ 27rZ, then 


lim 

6^0,R-^oo 





(fcizi, ...,Zn;(3) 


exists and is finite and real. Moreover, if f : C ^ C is a conformal map such that 
A = fizi), • • •, 4 = f{zn), then 


<Pc{z[,... ,z'^]f3) = Y[\f'{zj)\ sP “"^P(/c(^i,...,^n;/3)- 
i=i 
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Proof sketch. The beginning of the proof proceeds like that of Theorem 4.1 nntil 
eqnation (4.2), leading to the following equation: 




n 


exp 


b'=i 


5c{^i,...,^n},l<|S'|<n 


-XaR^S) ( 1 - cos ( /jfc 


\k&Is 


JJexp [-Aq; 5 ,k( 2 ;j )(1 - cos(3j)] , 
j=i 


where q:_r(S') = ficil ■ S G j,S'^ G diam( 7 ) < R), for S G {zi,..., Zn} with 
2 < |*S'| < n, and as,R{zj) = p.c (7 : ^ < diam( 7 ) < R, Zj G 7 , 2 ;a; ^ 7 Vfc 7 ^ j), and where 
Is denotes the set of indices such that /c G Is if and only if Zk G S. 

Note, in the equation above, the condition 151 < n in the hrst product on the 
right hand side; this condition comes from the fact that the term —XaR{S) with S = 
{zi,... ,Zn} is multiplied by 1 — cos{J2^=i Pk) = 0, where we have used the “charge 
conservation” condition \(3\ = J2j=i ^ 27rZ. 

For every j = 1,... ,n, using Lemma 8.1, we have that 


<^5,Ri^j) = hc(7 : ^ < diam( 7 ) < R, zj G 7 , Zfc ^ 7 V/c 7 ^ j) 
= ficil ■ 'kn > diam( 7 ) > 6, Zj G 7 ) 

+ /r ^(7 : m < diam( 7 ) < R, Zj e^,Zk ^ j) 

1 TTl 
= glogy + 


Now note that monotonicity and the fact that the Brownian loop soup is thin [17] im¬ 
ply that ajn,c{zj) = limi^^oo otm,R{zj) and ac{S) = lim^^oo ai?(5), for S G {zi, ...,Zn} 
with 2 < |5| < n, exist and are bounded. After letting R ^ 00 , the proof proceeds like 
that of Theorem 4.1, with D = D' = C. □ 


We have already seen the behavior of the 2-point function in the layering model in 
Section 3.4; the theorem below deals with the 3-point function. 

Theorem 4.5. Let zi,Z 2 ,Z 5 G C 6 e three distinct points, then we have that 


(j)c{zi,Z2,Z3-,/3i,P2,l33) = 


C3 


- Z2 


A,(/3i)+A7/32)-A,(/33) 






A,(/32)+AK/33)-A,(/3i) 2 


for some constant C3. 
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Proof. Theorem 4.4 implies that the 3-point function in the full plane transforms 
covariantly under conformal maps. This immediately implies the theorem following 
standard argument (see, e.g.^ [19])- We briefly sketch those arguments below for the 
reader’s convenience. 

Scale invariance, rotation invariance, and translation invariance immediately imply 
that there are constants Cahc such that 

^25 ^3) Aj 1^2;/^s) = 'y ^ CabcZl2 ^13^23 ’ (4-4) 

where Zij = \zi — Zj\ and the sum is over all triplets a,b,c > 0 satisfying a -|- 6 -|- 
c = 2(A/(A) + ^z(A) + ^z(A)) (the constraint on the exponents a,b,c follows from 
Theorem 4.4 applied to scale transformations). 

Now let / be a conformal transformation from C to C; / is then a Mobius transfor¬ 
mation and has the form f{z) = with f'{z) = ■ Letting = \ f'{zj)\~^, if 

5 = f{z), it is easy to check that Zij = Using this fact and Theorem 4.4, 

we have that 

0 c(AA 2 A 3 ; A, A, A) = 7^4/32 

_ [ Aim Aim Aii03)Ys^^ 

12 13 J /_^'^abc al2+bl2 a/2+c/2 f)/2+c/2 ' 

7l 72 73 

For this last expression to be of the correct form (4.4), the y’s need to cancel; this 
immediately leads to the relations a = 2 (A;(A) + ^/(A) ~ ^z(A)); b = 2 (A;(A) + 

AKA) - AKA)), c = 2(AKA) + AKA) - AKA)). □ 

5 Dimension of the exponential winding operator 

In this section we compute the conformal dimension of the exponential winding number 
operator for the BLS in the plane using non-rigorous path integral methods. 

The measure for a single Brownian loop of time length t rooted at x is given by 
the path integral for a free particle in two Euclidean dimensions, with the path y{T) 
constrained to begin and end at x: 

/ [cfy]e-yo<i^ \y\\ ( 5 . 1 ) 

Jy(0)=x 

where y = dy/dr, y,x are complex coordinates, and the path integral measure is 
Brownian normalized so that J^[(Py]e~^ /{2TTt). 

To construct the Brownian loop soup in the plane, we should integrate the root 
point X over the plane with uniform measure, integrate over t with measure dt/t, 
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and sum over n-loop sectors weighted by A" and divided by n\ since the loops are 
indistinguishable: 






n\ 

n=l k=l ' 


'52 


dtk 

tk 


(fxk 


i'yk{tk)=xk 


'ykio)=Xk 


[d'^yk] exp 



= exp 



(5.2) 


This result is closely related to (3.1) of [4], 

We wish to compute the 1-point function where (...) represents the av¬ 

erage over the distribution dehned by (5.2). Since is integer the result should be 
invariant under f3 ^ (3 + 27m. Furthermore for every conhguration there is a mir¬ 
ror conhguration where N —)■ —N, so is real and invariant under fd —)■ —fd. 

Because we are in the plane and charge conservation cannot be satished, the 1-point 
function will diverge as a power of the ratio of the cutoffs R/5, but the power tells us 
the dimension of the operator we wish to compute (c/. (3.3)). 

To compute we can insert in the path integral for each loop in 

(5.2). By translation invariance in the plane we can set z = 0. Then this insertion 
corresponds to adding a term i/d Jq dr 0(r)/27r to the single-particle “lagrangian” in 
the exponential of (5.1), where 0(r) mod 27r = arg(7/(r)) is the angular position of the 
path 7/(r) relative to the origin z = 0. Therefore 




exp 



ry{t)=x 


d^X 




'y(0)=x 


[d^y]e' 


Apart from the integral over t, the quantity in the exponential is related to the canonical 
partition function at inverse temperature f of a charged particle in the held of a single 
magnetic monopole of charge /9/27r at the origin (n-point functions involving products of 
exponentials would correspond to multiple magnetic monopoles with charge 

/dn/‘2n at the locations Zn)- To compute it, one can quantize the Hamiltonian for a 
particle in the monopole held, then take the trace in the energy basis [20], or perform 
the position-space path integral directly and obtain a result in terms of a sum over 
Bessel functions [21]. The hnal result is (see for instance [22]): 

/ l‘y{t)=x 

d^x / [d2j/]e-/o Liys/dr4>/2^ - 1 ) = -/9(27r - ld)/d>Ti^ 

Jy{0)=x ^ 2 


valid for 0 < fd < 27i and periodic in /? ~ 


= exp I -A 


Stt^ 


(d -1- 27m. From this we obtain 


(5.3) 
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in agreement with (3.3). 


6 Central charge and relation to the free boson 

The partition function (5.2) has a simple interpretation for the case A = 1/2. It 
corresponds precisely to the partition function for a free, massless, real bosonic held in 
two (Euclidean) dimensions: 

. ^ ^ rr ^ rdt^ _,rr 1 /■“dt /■ ,2 / , -/n, \ 

InZboson = --logdetid= - J y Tre ^ = - J — J d^x {x\e ^^\x) , 

where id = □ is the Laplace operator. Expressing the heat kernel using the 

path integral (5.1) completes the identihcation with (5.2). 

An analogous relation between the partition function of the discrete Gaussian free 
held (the lattice analog of the free held) and that of the random walk loop soup (the 
lattice analog of the BLS) can be derived easily; namely Z^Qpp = ^ 1/21 where 

ZpGFF die partition function of the discrete Gaussian free held in D with zero 
boundary condition, Z^^ is the partition function of the random walk loop soup in 
D with intensity 1/2, and \D\ is the number of vertices in D (see, for example [23], 
particularly Sections 2.1 and 2.2 and Exercise 2.5). This result may be related to work 
of Le Jan, who demonstrated that for A = /c/2, the occupation held for the BLS can 
be identihed with the sum of squares of k copies of a free held [16]. It suggests that 
for general A, the BLS might give meaning to the notion of a fractional power of a free 
held. 

Because the massless boson is a GET with central charge c = 1, the central charge 
of the BLS appears to be c = 2A for the case A = 1/2. Given the form of (5.2), one can 
express the partition function of the BLS for arbitrary A in terms of that for A = 1/2: 
Z\ = {Zii^) . The same relation can be rigorously proved to hold for the random 
walk loop soup mentioned earlier (see equation (2.2) in Section 2.1 of [23]). This leads 
to the conclusion that the relation c(A) = 2A between the intensity A and the central 
charge c of the BLS holds for all A. 

This agrees with the central charge of the SLE corresponding to the ensemble 
of BLS cluster boundaries. Indeed, conformal held theory considerations lead to the 
formula 

, , (3k — 8) (6 — k) 

^ Yk ’ 

but it is also known [2] that, when 8/3 < k < 4, the ensemble of cluster boundaries of 
a BLS with intensity 

(3k — 8)(6 — k) c(k) 

^ Ak ^ 
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is a CLE^. (For example, a BLS with A = 1/4 gives a CLE3, with central charge 1/2.) 

We note that most of the existing literature, including [2], contains an error in 
the correspondence between k and the loop soup intensity A. The error can be traced 
back to the choice of normalization of the (inhnite) Brownian loop measure which 
determines the constant in front of the log in Lemma 8.1 below.With the normaliza¬ 
tion used in this paper, which coincides with the one in the original dehnition of the 
Brownian loop soup [1], for a given 8/3 < k < 4, the corresponding value of the loop 
soup intensity A is half of that given in [2]. 

7 Conclusions 

There is still only a partial understanding of the rich connection between the conformal 
held theories studied by physicists and conformally invariant stochastic models such as 
the BLS or SLE. Physicists are often interested in CFTs dehned by a Lagrangian, and 
the correlation functions of primary operators, the spectrum of conformal dimensions, 
and the central charge are the main objects of study and interest. By contrast, confor¬ 
mal stochastic models are often dehned and studied using very diherent methods and 
with diherent goals. 

Here we have tried to take a few steps towards strengthening this connection. 
Using rigorous methods, we dehned a set of quantities in the BLS and proved that 
their expectation values behave like the correlation functions of primary operators. 
While we have not established it, we expect that these correlation functions may dehne 
a conformal held theory. If so, it has several novel features, such as a periodic spectrum 
of conformal dimensions. 

There are many basic questions that remain to be answered. Does this approach 
to the BLS in fact dehne a GET? If so, have we found the complete set of primary 
operators? What is the stress-energy tensor? Is the theory rehection positive and/or 
modular invariant? Is it unique in some sense? 

An even more ambitious set of questions relates to eternal inhation, albeit in 2-1-1 
dimensions. The model of [4] was proposed as a toy model for the late-time evolution 
of an eternally inhating spacetime, but the lack of analyticity in the 4-point function of 
the (analog of the) layering operator derailed it as a putative GET. Does the BLS solve 
this problem, and if so, could it help dehne a GET dual to eternal inhation? What 
type of object in 2-1-1 dimensional de Sitter space produces the BLS as its late-time 
distribution? Is there a natural generalization to higher dimensions? 

^We thank Greg Lawler for discussions on this topic. 
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8 Appendix. The Brownian loop measure: two lemmas 

In this appendix we prove two important lemmas which are used several times in the 
rest of the paper. The lemmas concern the /i^°'’^-measures of certain sets of loops, 
where, as in the rest of the paper, is the intensity measure used in the dehnition 
of the Brownian loop soup (see equation (1.1)). 

Lemma 8.1. Let z E C, then 

^ G 7 , (5 < diam( 7 ) < R) = ^ log 

5 0 

Proof. Since /i ^°'’^(7 : diam( 7 ) = i?) = 0, we have that 

p '°°^(7 : z E^,6 < diam( 7 ) < R) - /i '°°^(7 : G 7,7 ^ 7 C 

= p'°°P (7 : G 7 , diam( 7 ) > 5, 7 C B^^s) - : 2 ; G 7 , diam( 7 ) > i?, 7 C 5^,r), 

where B^^a is a disk of radius a around z, and 7 (/l B^^s indicates that the image of 7 
is not fully contained in B^^s- By the scale invariance of the last two terms are 
identical, so that 

/i '°°^(7 : z < diam( 7 ) < R) = : 2 : G 7,7 ^ B^^s, 7 C B^^r) 

= c\ogj (8.1) 

for some positive constant c < 00 , where the last equality follows from Proposition 3 
of [18] and the fact that the Brownian loop measure satisfies the conformal restriction 
property. In order to determine the constant c, we use the fact that, for any r > 0, 

^loop : z E ^,1 < diam(7) < e'’) = (7 : ^ G 7,1 < < e^’') . (8.2) 

(Equation (8.2) is essentially a consequence of Brownian scaling and can be proved 
using standard techniques. The interested reader can consult, for example. Appendix 
B of [23]). 

Following [18], we compute the right hand side of (8.2) using the dehnition of 
Brownian loop measure and translation invariance: 

p ^°°^’(7 : 0 G 7 ,1 < f., < e^") = 



c Jl 


27rR 


LzAb ■ 0 G -f})dtdA{z) 



[ ^ ^ e 7}) dt dA{z) 


2-Ri‘i ^ 0 '* \ jc dA{z) j dt 
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where denotes expectation with respect to a complex Brownian bridge of time 
length t started at the origin, and where, in the last equality, we have used the fact 


that 




because of scaling. 
[24], is 


The expected area of a “hlled-in” Brownian bridge, computed in 



so that 

; Z G 7,1 < = I- (8.3) 

5 

Using (8.3) and (8.2), we obtain 


: 2 ; G 7 ,1 < diam( 7 ) < e’’) = -. 

5 


Comparing this to (8.1) gives c = 1/5. □ 


Lemma 8.2. Let z E C and fc G Z \ {0}, then 


j iiQg winding number k around z,S < diam( 7 ) < R) 


1 

2n‘^k‘^ 


log 


R 

J' 


Proof. It is easy to check that the measure on loops surrounding the origin induced 
by but restricted to loops that wind k times around a given 2 ; G C, satishes the 
conformal restriction property. Therefore, as in the proof of Lemma 8.1, we have 


^ has winding number k around z,S < diam( 7 ) < R) 

= /i*"°^(7 : 7 has winding number k around <f. ^2,5, 7 C 

1 ^ 

= Cfelogy, 

for some positive constant Ck < 00 , where in the last equality we have used Proposition 
3 of [18]. In order to hnd the constants c^, we can proceed as in the proof of Lemma 8.1, 
using the fact that the expected area of a “hlled-in” Brownian loop winding k times 
around the origin was computed in [24] and is equal to l/27rfc^ for A; G Z \ {0} (and 
7r/30 for A; = 0). □ 
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